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Abstract 

Based on B — > D^K~K*° decay data, which supports factorization, we propose to extract the 
kaon axial form factor in the factorization framework. Experiment indicates that the K~K*° pair 
is produced by an axial current where only one out of three axial form factors is dominant. The 
axial form factor can be extracted by fitting the K~K*° mass spectrum with an ai(1260)-resonance 
plus QCD-motivated non-resonant contributions, which can be improved as data improves. 

PACS numbers: 13.25.Hw, 14.40.Nd 
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I. INTRODUCTION 



The three-body B -> DR-R*°, D*R~R*° and B -> DR-R°, D*R-R° decays were 
observed for the first time by the Belle Collaboration based on 29.4 fb _1 of data For the 
B -> D^R-R*° channels, one has a peak near threshold in the R R*° mass spectrum, 
which can be described by a dominant ai(1260) resonance contribution, although there is 
a possible peak at ~ 2 GeV. In addition, angular analysis finds R~R*° to be in the J p = 
1 + configuration, supporting the ai(1260) pole dominance picture. The fitted parameters, 
however, give larger a]" (1260) — ► R~R*° branching fraction |l], by a factor of 2 or more, 
than what is obtained from the three pion mass spectrum in r decays Q]. 

The peaking near threshold of the R~R*° mass spectrum suggests a quasi two-body 
process where the colinear R~R*° pair recoils against the meson. This suggests that 
factorization could be at work for such three-body decays. In our previous work on B° — > 
D(*) + K~K(*)° and B~ — > D^°K~K^° decays, two kinds of decay amplitudes arose due 



to different flavor structures that give rise to the D^*> meson [3|]. The B -> D^+R-R^ 
process involves only the matrix element (R~ R^°\V — A\0) , where the K~K^° is produced 
by a weak V — A current. For B~ -> D^°R~R^ , one has in addition a (K~K^°\V - 
A\B~) contribution, where B~ goes into K~K^° via a weak current. 

Under factorization, the B° — > D^ + K~K° decay amplitude is a product (Fig. E^a)) of 
the matrix elements (R-R°\V - A\0) and (D^ + \V - A\B ). From parity, K~ R° can only 
be produced by the vector current in (R~R°\V — A\0), which carries both J p = 0+ and 
1 _ components. The R~R° pair, however, should dominantly be in the 1~ state by isosp_in 
symmetry. This is consistent with experiment, which finds R~R° dominantly in 
By using isospin rotation, the kaon weak form factor (R~ R°\V\0) can be further related to 
the kaon electromagnetic (em) form factors in e + e~ annihilation, where much data exist. 
Without tuning parameters, we obtained the value for B(B — ► R~ K°) which is in 
good agreement with experiment 0. The predicted R~R° mass spectrum has a peak near 
threshold as a consequence of the kaon form factor, which can be checked by experiment. The 
prediction of R~R° to be in the 1~ state and good agreement with the observed branching 
fraction give strong support for factorization in B° — > R~ R° decay. 

As illustrated in Fig.d B -> D^ + R~R° decay bears some similarity with the baryonic 



yy isospin 

i-Gia. 



B — > D* pn decay, which was first measured by the CLEO Collaboration pj. Under 
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FIG. 1: (a) The K~K^ Q and (b) pn pairs produced by a current. The wavy line is the W boson. 

factorization, the B° — > D*~pn decay amplitude contains the matrix element (pn\V — 
A\0), which involves both vector and axial current form factors. We had also used isospin to 
relate the vector current form factor to the EM data. It was found that the vector current 
contribution could account for ~ 60% of the observed rate p. The predicted pn mass 
spectrum also exhibits near threshold enhancement. Unfortunately, this mass spectrum 
cannot be checked immediately. Furthermore, there is no data on axial current baryon 
form factors. We suggested instead that one could perform an inverse transform to obtain 
nucleon axial form factors once the B decay data becomes available. It would be interesting 
if fundamental quantities such as nucleon axial form factors can be extracted directly from 
B physics data. It should be noted that B° — > D^°pp spectra have been measured |2J and 
the ai-pole contribution could be important However, we cannot use them to obtain 
nucleon axial form factors, since these modes are dominated by transition contributions, i.e. 
A oc (pp\j\B°) (.D(*)°|j|0), in the factorization approach. 

Unlike the B° — > D*~pn case, we already have some data of the K~K*° mass spectrum 
for the B — > K~ K*° decay. In the factorization framework, the K~K*° pair in 5° — > 
D&+K-K* can be produced by both vector and axial currents. However, experiment found 
that K~K*° is dominantly in J p = 1 + configuration, implying a dominant axial current 
contribution. By neglecting contributions from the vector and the timelike component of 
the axial current (which would give 0~), it is possible to obtain the timelike K~-K*° axial 
form factors from B° — > K~ K*° decay data. 
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In this paper, our purpose is to demonstrate the extraction of axial form factor from 
data. Unfortunately, although more data is already available, at present the K~K*° mass 
spectrum is given only by combining D + K~K*°, D* + K~K*°, D°K~K*° and D* K~K*° 
modes. To stimulate further experimental studies, we generate mock data based on the 
existing one, from which an inverse transform is done to obtain the axial form factor. In 
the next section we formulate the factorization approach and define the K~-K*° axial form 
factors. The implications of angular analysis on the form factors will be examined, and 
the form factor is then parameterized for subsequent use. The fit results to mock data are 
obtained in Sec. IIH[ which is followed by discussion and conclusion in the last section. 



II. FORMULATION 



The relevant effective Hamiltonian is 

H cS = ^VebVZfaMOZfa) + c 2 (/i)0 2 c (/x)] , (1) 

where Ci(fj) are the Wilson coefficients and V c b and V u d are the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements. The operators Oi are products of V — A currents, i.e. 
0\ = (cb) V -A {du) V -A and O c 2 = (db) v ^ A (cu) v _ A , where (qq')v-A = 97 M (1 - 7 5 )<?'- 

We concentrate on the ~B° — > D + K~K*° decay mode, since under factorization it involves 
only a simple B° — ► D + transition and (K~K*°\V — A\0), giving the decay amplitude 

A(D+K-K*°) = ^y c6 U; d aa(i; + |(c6) y _ 4 iS Q )(ir-^* |(J W )^|0) , (2) 

where a\ = C\ + C2/N c with N c the effective number of colors if naive factorization is used. 
In practice, we may need to use the ai coefficient fitted from B — > Da\ decays. The matrix 
element (D + \(cb)v-A\B°) is the same as in two-body decay, and is parameterized by 

(B+(p D ) |U1 flV)) = (gr ~ ^f) (PB + PD)vF* D tf) + m| ~ ml q»F BD (q*) , (3) 

where q = ps — Pd — Pk + Pk*- We employ the Melikhov-Stech (MS) model |9| for the form 
factors F BD and Ff D . For (R- K*°\(du) v ^ A \0} , since K~K*° is produced by the V and A 
currents, in the K~K*° rest frame, the allowed angular momentum and parity configurations 
for K~K*° are J p = (U, 1 + and 1~, which transform in the same way as the axial vector 
current A^ (0~,1 + ) and the spatial part of the vector current V % respectively. A 
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convenient parametrization that manifests the possible J configurations is then given by 
i (K-(p K )K*°(j> K .,e K .) \{V - 0) = ie^e^p^ 2 * { ' r) 



VTIk + TTIk* 



[9,u-^)e^K + m K ,)A 1 (q 2 ) 



q 2 

9^u t- (Pk - Pk*) [e K * ■ q) ■ 

<T / m K + m K * 

+ ^ q,(e K , ■ q) A (q 2 ) , (4) 

where the form factors V, Aq and A\, A2 are induced, in the K~K*° rest frame, by the vector 
V 1 (1~), the timelike component of the axial vector A^ =0 (0~) and the spacelike component 
of the axial vector currents A 1 (1 + ), respectively. It should be noted that, once factorized, 
the angular momentum quantum number of the K~K*°, hence the allowed resonances, can 
only be J = 0, 1 to match the current. Any K~K*° or intermediate resonance component 
with J > 1 cannot be accommodated within factorization approach. 

Upon squaring the amplitude of Eq. 0, the decay rate would involve interference between 
the four form factors of (K~ K*°\V — A\0) . Simply knowing the K~K*° mass spectrum is far 
from being sufficient to constrain all the parameters. One thus needs further experimental 
inputs. Since we have parameterized (K~K*°\V — A\0) in accordance with the angular 
momentum configuration of the K~K*°, it is straightforward to see what information can 
be gained by studying the angular distributions. 

Denoting p£ the three-momentum of meson h in the K~K*° rest frame, and Ptin the 
B° rest frame, the helicity angle Qkk is defined as the angle between p* K , and — The 
angular distributions coming from V 1 (1~) and A^ =0 (0~) are 

oc sin 2 6 KK , (5) 



dYy . 2 



d cos 6kk 

oc constant , (6) 



dT Ao 



d cos 6 KK 

respectively, where we have used V and Aq to denote the corresponding contributions. 
Likewise, for A 1 (1 + ), we have 

dTAl oc l + bcos 2 6 KK , (7) 



d cos Qkk 

oc cos 2 Qkk , (8) 



dT A , 2 2 



d cos Qkk 

dT Al A* 
d cos Qkk 



oc cos 2 Qkk , (9) 
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FIG. 2: Fit to angular distribution by a(l + 6cos 2 Qkk) (solid), a constant c (dash), 
ism 2 9 KK (dash-dot) and ( cos 2 6 K k (dot), with (a, b, x 2 /n.d.f.) = (0.88,0.32,3.5/3), 
(c,x 2 / n -d-f.) = (0.99,4.0/4), (f, x 2 / n -d-f.) = (1.14,18.2/4) and (C,xV n -d-f-) = (2.0,39.3/4), re- 
spectively. The cos 2 and sin 2 Qkk distributions are clearly disfavored by data. 

from Ai, A 2 and their interference, respectively. 

As shown in Fig- El the data points scatter around the horizontal line <iN jd cos Qkk ~ 1- 
One concludes that the V (1~) and A 2 contributions are not favored by present data and 
can be safely dropped at this stage. On the other hand, angular analysis of subsequent 
K* — ► Kix decay suggests K~K*° is dominantly 1 + [l], therefore preferring A\ over A . 

It is interesting to understand how Ai, with dT a-l/ dcosOxK oc 1+6 cos 2 Qkk, can describe 
the data in Fig. EJ which seems quite consistent with a constant distribution. From 

dT tot ^ dT^i 
d cos Qkk d cos 6 KK 

oc J dM KK * \vr K A\VD\(F^ D ) 2 \A x \\^{\ + ^cm 2 9 KI ^ (10) 
oc a(l + 6 cos 2 Qkk) , (11) 

where = g 2 . We find that 6 in Eq. (fTT|) is small due to the factor in Eq. (fTUl) 

as a consequence of threshold peaking, resulting in a flattened parabola which mimics a 
constant. We reproduce the experimental data in Fig. El where one can see a clear threshold 
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FIG. 3: The observed W K~K*° mass spectrum combining D + K~K*°, D* + K~K* , D°K~K 
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and D*°K K*° modes. The experimental fit in 1] assumes only the ai(1260) resonance. The 
X 2 /n.d.f. is ~ 24/18 and the obtained nominal width of oi(1260) is 380 MeV. 

peak while large t = M\ K * is suppressed. The major part of the rate comes from the low 
Mkk* region, resulting in a suppressed in KK* frame. For higher M^k* where p|^* is 
large, the contribution to b is suppressed by the tail of the K~K*° mass spectrum. 

It should be emphasized that, although the present data disfavor any significant contri- 
butions from V, A 2 and A , it is possible that these contributions, however small, may show 
up in the future when one has sufficient amount of data. It is clear, however, that one can 
always use angular analysis to project out each form factor. As this is not yet the case, for 
this work, we concentrate only on the dominant Ai(q 2 ) term and neglect terms involving 
other form factors. The 5° — > D + K~K*° decay amplitude can now be simplified to 



A(D + K~K 



G) 



V cb V: d a^iq 2 ) A^q 2 ) (m K + m K .) 2p D ■ [e* K . - 



q)- (12) 



v/2 v q" 

As shown in Fig. El the data generally follows the experimental fit that assumes only the 
ai(1260) resonance, which we have now generalized into the A\{q 2 ) form factor. Note that 
the second peak at Mkk* ~ 2 GeV in the data has large statistical error, and is dubious for 
lack of known resonances. In any case, it would require more data to clarify. 

K~K*° can be produced via both resonant and non-resonant contributions. We pa- 
rameterize the A\ form factor by including a single ax (1260) resonance to account for 
the resonant contribution. For non-resonant contributions, we take cue from perturba- 
tive QCD (PQCD) |l0j . which states that at least one h ard gluon is needed to redistribute 
the large momentum transfer t in K~K*°. We then follow p] and expand in l/t n for the 
non-resonant part, with n starting at 1 to reflect this hard gluon exchange. For simplicity, 



and since data is still scarce at present, we take just two terms (n = 1, 2). We further 
impose a 180° difference between the phases of the complex expansion coefficients. This 
is to compensate for the (possibly) artificial rise at low t from the leading lit term. Such 
alternative signs are seen in the fits of nucleon EM form factors p . The non- resonant phases 
are taken as constant. The Ai form factor is then parameterized as 

n -1 



9s n^ai fai 



i<t> 



t 



X2 
t 2 



In 



t 



0.3 ; 



(13) 



t - m ax +im ai T(t) 

where m ai , f ai are respectively the mass and decay constant of ai(1260), g s the ai(1260) — > 
KK* coupling constant, and x\p and are the strengths and phase of the non-resonant 
terms. 

The di(1260) is broad and the width is not yet well known, and could range anywhere 
between 250 to 600 MeV It is not appropriate to treat the total width as constant, 

which is only valid in narrow width approximation. We use Q], 



T(t) 



ai 



P P n(m 2 ai ) ' 

where the constant T ai is the nominal width of ai(1260), and 



(14) 



PpiT 



1 + 



3 m 2 p 



(15) 



a, 



where, in addition to the phase space factor |p*| similar to what is used in [1|, we also take 
into account the amplitude squared of ai(1260) — > (p7r) s _ wave . The total width of ai(1260) 
is approximated by the ai(1260) — > (p7r) s _ wave partial width. 

The branching fraction of (2^(1260) — > K~K*° can be extracted by taking the ratio of 
the areas under the spectral functions of aj~(1260) — > all (approximated by a\ — > (p7r) s _ wave ) 
and a]"(1260) —> K~K*°, respectively: 



n 



r(t) 



p— 



n ( t _ m 2 i )2 + m 2 i r (t)i 



n 



m ai T a 



KK 



KK" 



n (t _ m 2 i)2 + m 2 ir(t) 2 

where the partial width of ai(1260) — > K~K*° is given by 

' 9s{m K + m K *) 2 p p7V (m 2 ai )] p KK *(t) 



(16) 
(17) 



KK 



•(*) 



57T 



rrin 



(18) 
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and 

pk 2 



Comparing to the definition of T(t), the constant inside the brackets of Eq. (fTHjl plays 
the same role as that of T ai in Eq. p4jl . i.e. the nominal value of the partial width of 
a^(1260) — > K~K*°. In the narrow width limit, T ai — > 0, we recover 



1 11 <-/ 9 

~ -<5(t - m 



(t-m2J 2 + m2 i r(t) 2 m ai r ai 



ai / ) 



such that 



IV « r(m^)5(t - r<) = r oi 5(t - m^) , (19) 
IW « r KK .(m 2 J8(t-m 2 J. (20) 

for r ai — > 0. Since in reality vtikk* > m ai the "branching fraction" i? of a^(1260) — > K~K*° 
is then defined by 

Our method for extracting the aj~(1260) — * K~K*° rate is similar to the one used by 
the CLEO Collaboration in Ref. Q]. In the fit to r~ — > z/ T 7r~7r 7r° which is dominated by 
t~ -> z/ r a^(1260), the a7(1260) -> if-if* "rate" (i? C LEO ~ 3%) was extracted by taking 
into account aj~(1260) — ► K~K*° and other processes that may contribute to the ai(1260) 
Breit-Wigner width. This is in contrast with the way the result was obtained in Ref. Q], 
where the ai(1260) was thought to account for the whole spectrum, and the corresponding 
factor i^Bdie is obtained by using ifeeiie = B(B° -> D + K-K*°)/B(B° -> L>+ar/(1260)) ~ 
15%. Only in the narrow width and m ai > tukk* limit that i?Beiie reduces to R defined in 
Eq. J2U). 



III. RESULTS 

Throughout this paper, we use a\ = 0.935, m ai = 1230 MeV, f ai = 229 MeV and the 
CKM matrix elements V ud = 0.975, V cb = 0.039 as in Ref. U. 

The mass spectra plotted in Fig. 01 is taken from Ref. |l(, which combines D + K~ K*° , 
D* + K~ K*° , D°K~K*° and D*°K~K*° modes. Our model parametrization in previous 
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section, however, is constructed only for D + R~R*°. Since our present purpose is only to 
demonstrate the feasibility of extracting the A\ form factor from B decay data, we shall 
make our own mock data on which to practice the extraction. 

Starting with Fig. El which is based on 29.4 fb _1 of data, we first transform from 
dN/dM KK * into dB/dM KK *, normalized by the measured B(D + R~R*°) = (8.8±l.l(stat.)± 
1.5(syst.)) x 10 -4 l|. Since the KEKB accelerator has already accumulated over 150 mil- 
lion BB pairs (~ 140 fb _1 ) by summer 2003, we take our mock data to be five times the 
original one (~ 29.4 fb _1 ), i.e. ~ 150 fb -1 , but keeping the central value for B{D + K~ R*°) 
unchanged. Next, we refine the resolution in Mkk* by rebinning two bins into three, again 
keeping the same central value and maintaining consistent statistical error. 

The second peak at M^k* ~ 2 GeV remains after rebinning, but no known resonances can 
account for it. We therefore remove it by hand. In order to do so, we first assume the data 
has a shape that roughly follows the curve given by some chosen set of parameters: g s = 1.4, 
r ai = 430 MeV, x x = 8 GeV 2 , x 2 = 1 GeV 4 , and = -40°, which give B(D + R-R*°) = 
8.4 x 10" 4 and x 2 /n.d.f. = 22.6/26. The branching fraction R of a^(1260) -»• R-R* given 
by Eq. (J2IJ is 4.3%, which is small compared to the value (8% - 15%) of Ref. |l|. With 
this curve as guide, we remove the second peak, but respect the statistical error in the 
corresponding Mkk* range. 
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FIG. 4: dB/dM KK * (in units of 1(T 4 ) mock data for B -» D + K~K*° assuming five times 
(~ 150 ftr 1 ) the original data |l( of 29.4 fb -1 . The solid curve corresponds to Eq. (|22|) . The dash 
curve is given by g s = 1.4, r ttl = 430 MeV, x 1(2) = 8(1) GeV 2 W and <f) = -40°. The dot-dash 
curve is given by g s = 0.89, r oi = 350 MeV, x 1(2) = 23.2(38.2) GeV 2 ^ and cf> = 201°. The vertical 
dot line indicates the r mass. 
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Having generated the mock data, we then redo the fitting. The best fit with x 2 /n.d.f.= 
19.2/26 is obtained with the following: 



g s = 2.41, r oi = 678 MeV, x x = 8.32 GeV 2 , x 2 = 3.07 GeV 4 , <p = -18.7°, (22) 



which give B(D + K~K*°) = 8.33 x 10~ 4 , and a larger branching fraction R = 8.18% due to a 
larger strong coupling g s . The fit result is shown in Fig.HJ The dominant contribution is from 
the non-resonant term, which gives ~ 54% of the rate. The ai(1260) resonance contributes 
only ~ 9% of the rate by itself, with the remaining ~ 36% arising from interference with 
the non-resonant term. The interference is constructive because the non-resonant phase 
(ft = —18.7° lies in the same quadrant as the resonance phase tan _1 [— m ai T(t)/(t — mlj], 
which varies between —90° to 0° for t > m\ . 

In Fig. 03 we show the A\ form factor with the parameters of Eq. (J22|) . The relative sign 
between the real and imaginary parts of the form factor is mainly due to the phase (f) = 
—18.7°, which is coherent with the ai(1260) resonance phase, and gives a negative imaginary 
part. The x\ term contributes ~ 131% to the non-resonant contribution and dominates A x , 
compared to only ~ 2% from x%, where we have built in destructive interference between x\ 
and xi parts. 

It is useful to investigate the x 2 behavior in the vicinity of the best fit values of Eq. (|22j) . 
We find that, within A% 2 = 1, the x 2 depends on g s and T ai roughly via the ratio g s /T ai , 
hence many "solutions" exist with differing a\ width. This is because the rate is dominated 
by the near-threshold contribution where the t — m 2 ai factor in Eq. (fT3*j) is small compared 
to m ai T(t), hence resulting in a resonance term that is roughly characterized by g s /T ai . 
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FIG. 5: Plot of |Ai| (solid), Re[Ax] (dash) and Im[Ai] (dot-dash) from Eq (JUJ). 
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Since this is inherent in our form factor model, one cannot determine g s and T ai separately 
even when the independent B° — > D + K~K*° spectrum has become available. One therefore 
needs independent input on ai(1260) resonance, which we discuss in the next section. 



IV. DISCUSSION AND CONCLUSION 



We see that independent input on ai(1260) width is needed to extract the A- A* axial 
A\ form factor. To clarify this point further, we note that the r — > u T K~K*° process 
is very similar to B° ->• D + K~K*°, with aiV cb {D + \(V - Ay\lf) in Eq. replaced by 
■u !/T 7 /t (l — 75)u T . Factorization in this case is basically exact, and A~A*° production comes 
only from the axial current. By us.ng the best fit valnes of Eq. Q, we obtain B(r - 
u T K-K*°) = 5.3 x 10~ 3 ~ 0.5%, similar to Ref. and is more than twice the measured 
rate of (2.1 ± 0.4) x 10" 3 Q. 

Our larger predicted r — > v T K~K*° rate can be understood as follows. Although both 
B(B° -> D + K~K*°) and B(t -»■ z/ r A~A*°) are proportional to |Ai| 2 , r -> v T K~K*° 
receives contributions only from the region of Mkk* < n^r- Our mock data implies ~ 60% 
of the & -> D+K-K* rate comes from M^k* < Tn T . The dominance of near-threshold 
contribution suggested by 5° — > D + K~K*° data translates to a large r — ► v T K~K*°. Thus, 
to reduce B(r — > z/ T A"^A"*°) from ~ 0.5% down to ~ 0.2%, one needs to reduce the relative 
contribution coming from near-threshold. We illustrate this by the dot-dash line in Fig. 
where only ~ 32% of the rate comes from Mkk* < m r and the peak is shifted towards 
higher t. This would have given the correct rates B(t — > v r K~K* Q ) ~ 2.1 x 10~ 3 and 
B(B° -> D + R-K*°) ~ 8.8 x 10~ 4 , but seems to disagree with the mock B° -> D + K~K*° 
spectrum. 

Recall, however, the comparison between the predicted spectrum of D + K~K° with that 
of D°K~K° in Ref. 0. Due to the additional B — > A" A" transition mechanism that is 
asymptotically 1/t 2 , the D°K~K° spectrum has a peak closer to the threshold (~ 250 MeV 
above threshold) than that of D + K~K° (~ 600 MeV above threshold). Since the mock 
data used in the present work is generated from published data that combines those of 
D + K-K*°, D* + K-K*°, D°K~K* and D* K~K*°, it is conceivable that the real data for 
B° — > D + K~K*° could be somewhat closer to the dot-dash line in Fig. which peaks at 
~ 400 MeV above threshold. In other words, the peak near the threshold of Fig. |H1 (or 
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Fig. 0J) could be due to the D^ K~ K*° modes, whereas the apparent peak at ~ 2 GeV in 
Fig. El (which we removed in FigEJ) might be some hint of the behavior of the D^*' + K~ K*° 
modes. The situation can be clarified by separating D + K~K*° and D°K~K*° with more 
data. We note further from Fig. El that \Af K (q 2 = (1.5 GeV) 2 )l ~ 1.3 is larger than 
\Ff K (q 2 = (1.5 GeV) 2 ) | ~ 0.7 as fitted from kaon EM data Ref. p. This gives further 
support of our conjecture that contamination of B — > K~ K*° in the mock data could 
be responsible for the enhancement of our r — > v T K~K*° rate. 

From the point of view of extracting the timelike K~K*° axial A± form factor, r — > 
v T K~ K*° is better because the production of KK* is purely proportional to A\. However, 
it has limited range in momentum transfer. In contrast, the 5° — > D + K~K*° mode allows 
us to probe higher t, but one has the associated factor aiFf D (see Eq. (II 2j) ). which means 
one is further subject to possible corrections to factorization. We have assumed certain 
values for aiFf D to extract A\. A more appropriate way to proceed would be to treat A\ 
and aiFf D (q 2 ) on equal footing, both as suitably parameterized functions to be determined 
from data. To be able to do so, one needs further independent data, such as an analysis of 
B — > Dp7r decay aimed at extracting B — > Da! (1260), which should be performed in a similar 
way to our B° — ► D + K~K*° study, including angular analysis and proper parametrization 
of iT-p axial form factor that respects PQCD for large M pn . There is, of course, the need for 
concurrent check on validity of factorization. This may then involve extraction of F^ D (q 2 ) 
from B° — > D + £~u. 

Our result reveals that the contribution from the ai(1260) resonance (~ 9%) is much 
smaller than that from the non-resonant part (~ 54%). While this underlines the importance 
of non-resonant part, the smallness of the ai(1260) resonance contribution is curious. Recall 
that in the B° — ► D + K~K° decay we have p-resonance contribution up to 40% [3j. It 
should be pointed out that, through SU(3) relations, the p resonance contribution in 

F KK 

is highly constrained by the well measured <fi resonance contribution in kaon EM data Q|. It 
is not clear if the smallness of the a\ resonance contribution is an artifact of having only one 
resonant term. On the other hand, the smallness of R ~ -Rcleo < 10% seems to support a 
suppressed ai(1260)-resonance contribution. It is interesting to see that, if we follow Ref. 
to take the ratio between the rate contributed from ai(1260) alone (9% of the total) and 
the rate of 5° — ► D + a^(1260) from Ref. 11|, we obtain R ~ 1.3%, which is about 9% 
of the Belle result -Rseiic ~ 15%. The discrepancy could be due to the breakdown of the 
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narrow width limit and the physical fact of itikk* < Wax, which are essential for -Reeiie to 
be identical with R defined in Eq. (|2ip. 

The case of B° — > D + p°7T~ is somewhat better in the sense that pn is dominated by the 
a\ resonance. However, there are many resonances that decay to the pit final state. One 
can extract B — ► Dai (1260) rate by suitably considering these resonant and non-resonant 
contributions. The branching fractions of B(B° -> D+a^(1260)) = (0.60±0.22±0.24)% and 
B(B~ -> L>%(1260)) = (0.45 ± 0.19 ± 0.31)% were obtained by an analysis of B -> Dtxtxtx 
decays based on ~ 200 pb _1 of data The large systematic error is dominated by uncer- 
tainties in fitting and are estimated by considering alternate backgrounds. The mass and the 
width of ai(1260) were taken as a priori in the extraction, which may no longer be appropri- 
ate when the data has improved by more than three orders of magnitude by the B factories. 
To reduce the effect caused by uncertain properties such as the width and sub-decay modes 
of ai(1260), an amplitude that takes into account as complete substructures of ai(1260) as 
possible is needed. Furthermore, as we have seen from the fit to the B -> D + K~ K*° data, 
a non-resonant part that respects PQCD at higher t may contribute significantly to the re- 
sults, which must also be taken into account to complete the construction of the amplitude. 
Of course, the cost would be to boost the number of fit parameters, but fortunately one now 
has enormous amount of data. 

As noted in the end of the previous section, the ai(1260) contribution depends on g s 
and r oi roughly through the ratio g s /T ai . Thus, B° — > D + K~K*° mode is not so powerful 
for determining the resonance parameters. One needs independent determination of these 
parameters, especially r oi from other sources. Therefore, it appears that, to extract the 
axial K-K* form factor Ai and understand properly the ax (1260) component, one needs to 
perform a combined analysis of B° — ► D + K~ K*°, D + p°n~, D + t~v, and r — > v T K~K*°. 

In conclusion, based on what is revealed by 29.4 fb" 1 data from Belle, we have illustrated 
how to extract the timelike K~-K*° axial form factor Ai(M|- K ») from B° — > D + K~K*° data. 
The method can be easily generalized to other form factors by angular analysis when more 
data becomes available. With > 300 fb _1 data already accumulated by KEKB and PEP-II, if 
the case for factorization in such three-body B decays is further strengthened, the B Factories 
may be promising for the study of hadronic form factors that may be otherwise inaccessible. 
Hadronic form factors have played an instrumental role in the formulation of nuclear and 
particle physics, and further insight may perhaps be gained from quantities such as the 
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K-K* axial form factor. Our study already indicates that the ron-resonant contribution, 
needed to account for large t behavior, is rather significant. The current discussion also 
illustrates how nucleon axial form factors can be extracted from B° — > D^~pn decay once 
the spectrum becomes available. 
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